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ABSTRACT 

The helical magnetorotational instability of the magnetized Taylor- Couette 
flow is studied numerically in a finite cylinder. A distant upstream insulating 
boundary is shown to stabilize the convective instability entirely while reducing 
the growth rate of the absolute instability. The reduction is less severe with larger 
height. After modeling the boundary conditions properly, the wave patterns 
observed in the experiment turn out to be a noise-sustained convective instability. 
After the source of the noise resulted from unstable Ekman and Stewartson layers 
is switched off, a slowly-decaying inertial oscillation is observed in the simulation. 
We reach the conclusion that the experiments completed to date have not yet 
reached the regime of absolute instability. 



Subject headings: accretion, accretion disk- 
MHD — methods: numerical 



-inst ability — (magnet ohydro dynamics : ) 



1. Introduction 



The magnetorotational instability (MRI) is prob ably the main source of turbulence and 
accretion in sufficiently ionized astrophysical disks (jBalbus &: Hawleyl Il998l ). Due to this 
crucial role in astro physics, substantial efforts have been spent worldwide to observe MRI in a 
laboratorv setting Jji et all boOll: ICoodman fcjilbood : iNoguchi et aPbooi Isisan et al.lbood : 
Velikhov et aDbooef L but MRI has never been conclusively demonstrated in the laboratory. 
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Most experiments have been done in cylindrical geometry with a background flow that 
approximates the ideal Couette rotating proflle: 



Q = a + b/r"^ 



(1) 



where a = (f22r| — i7ir^)/(r2 — r^) and b = rjr|(r2i— i72)/('"2~''"i)' ^"^^ ^2 are the rotation 
speed of the inner and outer cylinder and ri and r2 are the radius of the inner and outer 
cylinder, respectively (see Fig. [1]). For axially periodic or inflnite magnetized Taylor- Couette 
flow, MRI-like modes have been shown theoretically to grow at much reduced magnetic 
Reynolds number Rem = fii?"i(^2 — 1^1) /v and Lundquist number S = Va,o'^i('^2 — fil/v in 
the presence of a combination of axial and current-free toroidal fleld 



(2) 



than the standard MRI (SMRI) with purely axial magnetic fleld (IHollerbach fc Riidiger 



2005 



Riidiger et al.l 120051 ) . Here the cylindrical coordinates (r, y9, are used. and (3 
are constants. The Alfven speed is deflned as Vaa = B^/ y/Aiip. rj and p are the magnetic 
diffusivity and density of the fluid, respectively (see Fig. [T]). 

The Potsdam ROssendorf Magnetic Instability Experiment (PRO MISE) group claimed 



to have obseryed this kind of helical MRI (HMRI) experimentally (IStefani et al.l 12006 



Riidiger et al.l l2006l : IStefani et al.l 120071 ) . However we have shown that the wave pattern 
observed in PROMISE is not a global instability, but rather a transient disturbance some- 
how excited by the Ekman circulation and then transiently amplifled as it propagates along 
the background axial Poynting flux with nonzero group and phase velocities, but is then 
absorbed once it reaches the jet formed at midheight between two neighboring Ekman cells 
(ILiu et al.ll2007l ). PROMISE group have accordingly updated the experimental facility to 
PROMISE II to allow for two split rings at both endcaps: the inner ring attached to the 
inner cylinder and outer ring attached to the outer cylinder. If the width of the inner ring is 
chosen appropriately ~ 0.4(r2 — ri), the magnetized Ekman circulation could be signiflcantly 
reduced, therefo re removing one of the possible disturbance sources, i.e. the unsteady jet 
(|Szklarskill2o"o7l ). 



As with other examp les in the literatires, such as drifting dynamo waves (ITobias et al. 



19981 : iProctor et al.ll2000l ). it is of vital importance to distinguish absolute instability from 
convective instability in a traveling wave experiment like PRO MISE. It is also an essen- 
tial ingredient of the threshold prediction for the Riga dynamo (IGailitis et al.ll2008l ). For 
a traveling wave the positivity of the growth rate implies only an ampliflcation of the per- 
turbation as it moves downstream. In one case, despite the movement of the wave packet, 
the perturbation increases without limit in the course of time at any point flxed in space; 
this kind of instability with respect to any infinitesimal perturbations will be called absolute 
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instability. In the other case, the packet is carried away so swiftly that at any point fixed in 
space the perturbation te nds to zero as t ^ cxd; this kind will be called convective instability 
( iLandau &: Lifshitzl 119871 ) (see the details of For PROMISE II, it appears that under 
the experimental conditions the second kind occurs. A recent preprint also highlights the 
importance of the distinction betwe en absolute and convective instabilities in the context of 
the HMRI JPriede fc Gerbethlboosh . 



In a Taylor- Couette experiment bounded by insulating endcaps, iTobias et al.l (Il998l ) 
have pointed out that without any external disturbances except a small initial disturbance 
needed as a seed for the instability, the distant upstream insulating boundary acts as an 
"absorbing" boundary while the characteristics of the downstream endcap is unimportant. 
Due to this absorption the convective unstable state ca nnot be sustained b y a uniform driving 
force, therefore this unstable mode eventually decays fiTobias et al.lll998l ). This driving force 
is not the noise mentioned before, but the power to drive the instability, which in the usual 
Taylor- Couette experiments can be quantified by the magnetic Reynolds number Rem- This 



conclus ion has been rigorously demonstrated in th e very resistive lim it in §11. C of iLiu et al. 



(l2006al ) using a perturbative approach and §11. D of iLiu et al.l (l2006al ) using a modifed WKB 
analysis, showing that the insulating endcap entirely stabilizes the HMRI mode, which is a 
convective unstable mode given the parameters of the PROMISE experiment. 

The absorbing boundary is essential to the development, regardless of how distant it may 
be. The larger height onl y defers the time w hen we have to wait for the boundary-induced 
dissipation to dominate ( Tobias et al.l 1998). On the other hand, if Re^ exceeds a higher 
threshold -Rcm.f, the dri ving force of the s ystem overcome the dissipation and a globally 
unstable mode appears (ITobias et al.lll998l ). Therefore in a bounded system the unstable 
mode appears at Rera,i rather than Rem^c, where Rera,c is the critical magnetic Reynolds 
number for the ons et of the convective unstable mode without the "absorbing" boundary. 
Tobias et al.l (119981 ) has showed that in the presence of an "absorbing" boundary and large 
h, a global unstable mode appears when 

Re^ ^ i?em,f = i?em,a + 0{h~^) , 

where Rem,a. is the critical magnetic Reynolds number corresponding to the onset of the 
absolute instability without the "absorbing" boundary. 

This has raised a big obstacle for people to observe absolutely unstable HMRI in the 
laboratory. The advantage of HMRI itself, i.e., unstable with a low critical Reynolds number 
(3 orders lower than the SMRI) confiicts with the necessarily high threshold of the onset of an 
absolute HMRI mode, i.e., excited at a reasonably high critical magnetic Reynolds number, 
thus high Reynolds number, which would result in much more severe end-effects than people 
had expected. Moreover the fact that the critical Lundquist number must usually increase 
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together with the magnetic Reynolds number and high ratio of toroidal-to-poloidal magnetic 
field requirement {j3 > 1) would even worsen the situation. 

We also find that by nonlinear numerical simulation the insulating endcap reduces the 
growth rate of the absolute instability somewhat. The higher the height h is, the less the 
growth rate is reduced (Tabled]). 

In a typical experiment, the experiment is, however, highly likely affected by small exter- 
nal noise either from a physical cause or experimental imperfection such as the misalignment 
of the cylinders. If the system is convectively unstable, i.e., disturbances grow as they move 
downstream, noise would sustain s tructures in the system even if no global mode is unstable 
(jDeisslerlll987l : iProctor et al.ll200Cll ). In the present paper, we show by numerical simulations 
that the perturbations from the unstable magnetized residual Ekman layer and Stewartson 
layer at the upper endcap would play the role of "noise" generator, though this perturba- 
tion level is reduced with increasing axial magnetic field (jLiull2008al ). What is observed in 
PROMISE II turns out to be a noise-sustained convective traveling wave, not the absolute 
unstable mode. 



This paper is organized as follows: [|2] presents the wave packet analysis in a unbounded 
cylinder, which is the basis of the following sections. We report the nonlinear simulation 
results with partially conducting boundary conditions of PROMISE II experiment in ^ The 
final conclusions and implications to the HMRI experiments are given in §H 



2. Wave Packet Analysis in an Unbounded Cylinder 

Assuming a cylinder of infinite height h, kz is a continuous variable. Let the gap width 
be fixed and finite, so /cj. = 7r/(r2 — ri). We define the total wavenumber K = ^/k"^ + k1 and 
the growth rate 7. 

Since the fast growing mode is the dominant mode, here we focus on waves with vertical 
wavenumber kz close to that of the fastest growing mode, The range of values of kz lies 
near the point for which 7(^2) is a maximum, i.e. d'j/dkz = ai kz = k^ [as seen from Fig. [2] 
(a)]. Let a slight perturbation occurs near the middle of the fiow {z ~ 0) in the format of a 
wave packet as follows: 

Briz, t = 0) = 60 exp (^-Jj2^ exp(ifc°z) , (3) 

where we have used the envelope exp(— 2;^/2L^) to confine the perturbation around the 
central part of the cylinder, where L ~ 0{h). In the course of time, the components for 
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which 7(^2) > will be amplified, while the remainder will be damped. The amplified 
wave packet thus formed will also be carried downstream with a velocity equal to the group 
velocity du/dk^ of the packet, where uj = Mcj + i'j and M.uj is the real part of the frequency; 
since we are now considering waves whose wave numbers lies in a small range near the point 
where d'y/dkz = 0, the quantity 



Vg = du/dK = d{Ruj)dkz 



(4) 



is real, and is therefore the actual velocity of propagation of the packet. This downstream 
displacement of the perturbations is very important, and causes the complications of absolute 
instability v.s. convective instability. 



We can approximate the dispersion relation like (Fig. [2]): 

k 

Rio = muj{kz) = ; 



7 = l{kz) = 1^ 



k 



0n2 



(5) 
(6) 



in which 



1 d 



^^{r^Q)'^ = A{l + Ro)Q^ and Ro = l/2dlnQ/d\nr = a/Q-l is the Rossby 
number. We know 7 = when kz = 0. Thus a = 27*^//c°^. And in order to simplify the 
derivation, we assume K ^ constant from now on (though this is not a good approximation, 
we can get some insightful results from this simple approximation). From Eq. [5l we get 



- k/K. 
At later time t > 



Br{kz, t) 



= Br{kz, 0) exp{-f{kz)t + iMuj{kz)t) 
5;(A;„0)exp{[70 - - A;°)2]t + m^t} , 



(7) 



if we define D = + at, the result can be expressed as: 



Br{z, t) = bo— exp(7°t) exp 



[z + v,ty 

2D^ 



exp[zA;°(z + Vgt)] . 



In Eq. [8], As t — * 0, Eq. [8] can be simplified as: 

Br{z,t) = 6oexp(7°t)exp[2A;°(2 + 1/gt)] 
which is a "transient" growing phase. As t — > 00, 

L 



Br{z,t) = bo exp 



7 



9 
2a 



exp[iA;°(^ + Vgt)] 



(9) 



(10) 



Obviously, If 7" < 7a = Vg /'^o'^ we will get convective instability, that is, it starts with 
a transiently growing phase (Eq. [9]), followed by a phase asymptotically decaying to zero 
(Eq. [To]). If 7" > 7a, we will get absolute instability. 
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3. Noise-Sustained Convective Instability in PROMISE II Experiment 



In order to reduce the undesirab le effects induced by the endcaps and also the accom- 
panying hydromagnetic asymmetries, ISzklarskil (120071 ) have proposed to split both endcaps 
into two rings which are attached to both cylinders and found that if the width of the inner 
ring is chosen to be OAD (see Fig. [T]), where D = r2 — ri is the gap between the inner 
and outer cylinder, the magnetic energy in term of 6^, where is the perturbed azimuthal 
magnetic field, is minimized. Therefore the magnetized Ekman circulation is significantly 
reduced, leading to a satisfactory ideal Couette state (Eq. [1]) in the bulk flow. PROMISE 
has been accordingly updated to PROMISE II adopting this idea. 



While we have confirmed their conclusions (Fig. [3]) (Please note that in lSzklarskil (120071 ). 
this conclusion is derived with j3 = 0, i.e., no background toroidal magnetic field, while 
our simulation results show that this conclusion is also v alid with nonzero /?), here we 
report nonlinear simulations with the ZEUS-MP 2.0 code (iHayes et al.ll2006l ). which is a 
time-explicit, compressible, astrophysical ideal MHD parallel 3D code, to which we have 
added viscosity, resistivity (with subcycling to reduce the cost of the induction equation), 
and partially conducting boundary conditions (ILiu et al.l 120071 ). for axisymmetric flows in 
cylindrical coordinates (r, ip, z). It has been demonstrated that the finite conductivity (r^cu = 
1.335 X lO^cm^s^^) and thickness of the copper vessel are important, and this noticeably 
improves a greement with t he measurements compared to previous much simplified boundary 
condition (ILiu et al.l 120071 ). Please note that in this paper = VL2/VL1 = 0.26, rather than 
fi = 0.27 rep orted in previ o us wo rk. The parameters of PROMISE II as reported in or 
inferred from IStefani et al.l (120081 ) are used: gallium density p = 6.35 gcm~^, magnetic 
diffusivity rj = 2.43 x 10^ cm^ s~^, magnetic Prandtl number Pr^ = u/rj = 1.40 x 10~^; 
Reynolds number Re = Qiri{r2 — ri)/v = 1775; axial current = 6000 A; toroidal-coil 
currents = 0, 50, 75, 120 A; and dimensions as in Fig. [H 

For comparison, we start with purely hydrodynamic (unmagnetized) simulations (Fig.Hj). 
From Fig. H] (a), after splitting the endcaps into two rings, the two big Ekman cells are di- 
vided into four small er cells and loca lized near the endcaps. Compared to the simulation 
results of PROMISE (ILiu et al.ll2007l ). there is not an flapping "jet" near the mid-plane as 
in the usual Ekman circulations. This removes the possible noise from this unsteadiness. 
However from Fig. HI (b), there are some perturbations near both endcaps, which supply 
the possible sources of noise in the sys tem. Thes e perturbations are resulted from unsta- 
ble Ekman layer and Stewartson layer (lLiull2008al ). The magnitude of this noise is around 
±0.2 mm s~^. As we will see later (Fig. [5]), this unsteadiness is reduced by increasing axial 
magnetic field jGilmanlllQTll : iLiulboOSal l. 



Figure. E] displays vertical velocities near the outer cylinder in simulations corresponding 
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to the experimental runs of IStefani et al.l (120081 ) for several values o 



I^p. A wave pattern very similar to that in the experimental data (IStefani et al.l l2008l ) is 



the toroidal current, 



seen. Since now there is no jet, the traveling wave is propagating to the bottom endcap and 
abs orbed there whi le in the old PROMISE experiment, the traveling wave disappears at the 
jet (iLiu et al.l 120071 ). We also notice that the perturbation near the upper endcap weakens 
with strong axial magnetic field. This coul d be expla ined by a more stable magnetized 
residual Ekman layer and Stewartson layer (lLiull2008al ) . Both the weakening of the noise 
sources and disappearance of the amplifying mechanism leads to a rather steady state with 
= 120 A. 

It is highly possible that there is much noise in the real experiment due to some experi- 
mental imperfection such as misalignment and in the numerical simulation such as numerical 
noise. Also the noise could result from physical causes such as the unsteady Ekman layer or 
Stewartson l ayer. These no i ses w ould cause a noise-sustained convective instability in the 
system as in iProctor et al.l (120001 ) . The continuous impulse from the noise sources would 
have the system always in the state of "transiently growing" phase (Eq. [9]). This results 
in similar wave patterns as the ones from the primary instability wit hout noise, whi ch are 
observed in PROMISE and PROMISE II experiments and simulations (ILiu et al.ll2007l ). The 
noise-induced w ave pattern is always susceptible to noise-ind uced disruption a s discussed by 
Deisslerl ( 119871 ). That is exactly wh at we found here and in iLiu et al.l (120071 ). We can see 



this point more clearly by following iLiu et al.l (120071 ): performing a simulation that begins 
with the experimental boundary conditions until the traveling waves are well established, 
and then switches abruptly to ideal-Couette endcaps (Fig. [H]). After the switch, the traveling 
waves disappear after one axial propagation time and slowly decayin g inertial oscilla tions 
(asymptotically to zero) result. The main difference in results between iLiu et al.l (120071 ) and 
the present simulation are: (1) there is no jet, thus the traveling waves are absorbed near 
the bottom endcap both before and after the switch; (2) there is no change of wave speed 
associated with the switch since the background state does not change much before and after 
the switch. We reach the conclusion that even after the endcaps are split into two rings as 
in PROMISE II, the wave patterns observed in the experiment are not global instability, 
but rather noised-sustained convective instability. The similarity between these "inertial- 
oscillation-induced wave" after the switch and the earlier noise-sustained "MRI waves" in 
the simulation or "MRI-type waves" observed in the various versions of PROMISE stems 
from the physica l nature of HMRI that HMRI is a weakly destabilized inertial oscillation 



(iLiu et al. 



2006al ). More importantly, this similarity supports our conclusion in another as- 



pect: the frequency and wave number selection mechanism for a noise-sustained structure is 
determined by a linear mechanism, thus resembling the properties of the primary instability. 
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4. Discussion 

In this paper, nonlinear simulations of the helical magnetorotational instability in a 
magnetized Taylor- Couette flow are performed. The geometry mimics PROMISE II ex- 
periment wit h endcaps split into two rings. The partially conducting boundary condition 



introduced in lLiu et al.l (120071 ) is used. The waves patters change with applied magnetic field 
as in the experiment. However via numerical tests, we find that the wave patterns observed 
in PROMISE II experiment are not due to a global instability, but rather a noise-sustained 
convective instability. 

The importance of the distinction between absolute and convective instability in a 
bounded system with broken reflection symmetry is discussed. The addition of the toroidal 
magnetic field breaks the axial symmetry of the system. In such cases, the effects of distant 
upstream insulating boundaries on the absolute instability differs remarkably from the ones 
on the convective instability. The insulating endcap would only reduce the growth rate of the 
absolute instability, but would stabilize the convective instability entirely, however distant 
it may be. For the absolute instability, the more distant insulating endcap would less reduce 
the growth rate, while for the convective instability the more distant endcap would only have 
the system wait longer for the dissipation due to the "absorption" boundary to dominate. 
These discoveries cast great obstacles for people to observe the helical magnetorotational 
instability in the laboratory: An absolute HMRI is needed to observe the global unstable 
mode in the experiment. 

Unfortunately it is not easy to derive the critical magnetic Reynolds number -Rcm^f of 
the absolute HMRI analytically in a bounded system. However we can get a rough estimate 
of i?em a; ^-C, the critical magnetic Reynolds number of the absolute HMRI in an unbounded 
system, by wave packet analysis and the approximate dispersion relation from Fig. [2l 
From Fig. [2l we derive the group velocity Vg ~ 1.08 cm s~^, 7° ~ 0.31 s~^, ~ 0.52 cm~^ 
and a ~ 2.29cm^s~^. Therefore 7*^ — Vg/2a ~ 0.05 s~^ > 0, which corresponds to an 
absolute HMRI instability with Rem,ii ~ 0.07. We therefore conjecture that Rera,i = -Rem,a + 
0{h~^) > 0.07 in PROMISE II. The critical magnetic Reynolds number is somehow one order 
of magnitude lower than the standard MRI, but still requires Reynolds number Re ~ 10^. 
Therefore we need to rotate the cylinder typically with more than one hundred rpm. Such 
rotation rates are of course achievable, however with such a Reynolds number the advantage 
of HMRI with much lower Reynolds number, thus much lower end-effects, is not so great 
as people had expected. Moreover in most HMRI unstable modes /? > 1 is preferred, this 
suggests a toroidal magnetic field typically at ~ 1,000G, which requires axial currents 
> 10^ A inside the inner cylinder. This is a big engineering challenge in itself. The technical 
constraints prevent us to try to find the threshold by nonlinear numerical simulations such 
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as: (1) the current code can not afford large Reynolds number (~ 10^), which is required 
for HMRI to enter the absolutely unstable regime; (2) from the global linear calculation, the 
HMRI mode is stabilized if an artificially low Reynolds number like ~ 10^, which could be 
afforded by the current code, is employed. 



The non-axisymmet r ic m = 1 modes a re observed in the experiments (jStefani et al. 



20061 : iRiidiger et al.ll2006l : IStefani et al.ll2007l ) . Unfortunately since the simulations presented 
in this paper are all axisymmtric, the possibility to study this important mode is excluded. 
The extension of the current work to 3D will be the subject of the future study. Ruediger 
et al have already done some excellent work on this issue and found that given PROMISE 
parameters the n onaxisymmetric HMRI modes are always ha rder to be excited than the 
symmetric mode (IRiidiger et al.ll2005l : IRiidiger fc Schulta 120081 ). 
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20.3 


40.6 


81.2 


periodic 


Growth Rate 73"-'^ 


0.27 


0.58 


0.82 


1.06 



Table 1: Influence of the height h upon the growth rate 7 of the absolute instability in a 
bounded cylinder, ri = 7.1cm, r2 = 20.3 cm, VLi = 400 rpm, ^2 = 53.3 rpm, = 500 G, 
Belri) = 1 kG, the height h = 27.9 cm, 55.8 cm and 111.6 cm; the material properties are 
based on gallium: rj ^ 2000cm^s~^ and p ~ 6gcm~^, which give Re^ = 2 and S = 2.7, 
no explicit viscosity present. Th e simulations are perform e d using a modifi e d version o f 



no expncit viscosity present. 1 n e simulations are periorm e a using a moain e a version 1 
the astrophysical code ZEUS2D jstone fc NormanI ll992al E]: ILui et all boOfiljjLh] boOSbl ). 
The boundary conditions adopt the one introduced in §11. D of iLiu et al.l (l2006a[ ). Please 
note that no-slip boundary conditions are employed on all applicable boundaries and ideal 
Couette state (Eq. [T]) is enforced at both endcaps in order to remove the Ekman circulation 
and possible disturbances induced by this boundary layer effect. The one labeled "periodic" 
uses vertically periodic boundary conditions with periodicity length h = 27.9 cm. 
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IV 



III 



III 



II 



IV 



IV 



IV 



01 



9 



8.01 cm 
4.32 cm 

"r 

■~^4.11 cm 



16 cm 



Fig. 1. — Computational domain for simulations of PROMISE II experiment. Region (I): 
Inner copper cylinder, angular velocity Qi. (II): outer copper cylinder, (HI): liquid 
gallium; (IV): vacuum. Thick dashed line: insulating inner ring, corotating with the inner 
cylinder. Thick dash-dot line: insulating outer ring, corotating with the outer cylinder. The 
junction of these two rin gs lies at 40% of the gap {D = r2 — ri) between the inner and outer 
cylinder (jSzklarskil 120071 ) . Dimensions: ri = 4.0 cm; r2 = 8.0 cm; h = 40.0 cm; dwi = 1.0 cm; 
dwii = 1.5cm; fii/27r = 3.6 rpm; fl'?, /27i = 0.936 rpm. Note that u = p,y,/Q-\ = 0.26, rather 

' 1 " 



than At = 0.27 used in previous work (iStefani et al.ll2006l : iRiidiger et al.ll2006l : iLiu et al.ll2007 



Stefani et al.l 120071 ). The exact configuration of the toroidal coils being unavailable to us, 
six coils (black rectangles) with dimensions as shown were used, with 67 turns in the two 
coils nearest the midplane and 72 in the rest. Currents were adjusted to reproduce the 
reported Hartmann numbers Ha = B'^ri/ y/p/J^fjU. 
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Fig. 2. — (a) Growth Rate 7; (b) Real Frequency Mcj. * Linear Calculation, - Approximation 
by Eq. [5] and Eq. El ri = 4.0 cm, r2 = 8.0 cm, Qi = 101.25 rpm, Q2 = 26.325 rpm, Bz = 
220.5 G, (3 = 4.0; the material properties are based on gallium: t] = 2.43 x lO'^cm^s^^, 
u = 3.4 X 10~'^cm^s~ ^ and p = 6.35 g cm~'^. The calculations are performed using a code 
(IGoodman fc Jill2002l ) adapted to allow for a helical field. Vertical periodicity is assumed, 
but the radial equations are solved directly by finite differences with perfectly conducting 
boundary conditions (§11. B of iLiu et al.l (j2006al )). 
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Fig. 3. — Azimuthal velocity v.s. Radius r. Re = 1775, /3 = 3.81 and = 75 A. Solid 
line, ideal Couette state; +, 1.31 cm; *, 2.72 cm; □, 13.95 cm. 
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Fig. 4. — Purely hydro dynamic (unmagnetized) simulations. Left: Time-averaged poloidal 
flow stream function \E'; Right: (color) Axial velocities [mms~^] versus time and depth 
sampled at r = 6.5 cm, for the parameters of the PROMISE II experiment without any 
magnetic field. Note height increases upward from the bottom endcap. No-slip velocity 
boundary conditions are imposed at the rigidly rotating endcaps. The steady part of the 
resulting Ekman circulation is suppressed in right panel by subtracting the time average at 
each height. 
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Fig. 5. — (color). Axial velocities [mm s~^] versus time and depth sampled at r = 6.5 cm, 
for the parameters of the PROMISE II experiment with toroidal currents as marked. No- 
slip velocity boundary conditions are imposed at the rigidly rotating endcaps. The steady 
part of the resulting Ekman circulation is suppressed in these plots by subtracting the time 
average at each height. The waves appear to be absorbed near the bottom endcap. 
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Fig. 6. — (color). An extended version of the case I^ = 75A shown in Fig. [5] but without 
subtraction of the time average. After t = 360 s, the no-shp boundary condition at both 
endcaps is switched to an ideal Couette profile (Eq. [1]). A slowly decayed inertial oscillation 
is resulted. 



